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SUMMARY 

Two mechanisms of drag reduction for flow over flat plates are 
investigated. The first mechanism employs Bushnell’s hypothesis that 
compliant walls produce drag reduction by interfering with the formation 
and bursting of the turbulent spots in a turbulent boundary layer. It 
is shown that the required amplitudes and frequencies of the wall motion 
might be achieved by using slightly curved walls and the resulting large 
amplitude motions of snap buckling of the walls. A simple structural 
model of an arch is used in the analysis. It is shown from the exact 
solution of the arch problem that the required structural motions can be 
obtained for compliant walls constructed from materials like mylar. In 
obtaining the solution, a new asymptotic method, which we call the method 
of rational functions, was developed for analyzing jump phenomena. In 
the second mechanism the delay of transition from laminar to turbulent 
flow by driven walls is studied. It is shown for Poiseuille channel flow 
whose walls are driven by a periodic traveling wave that a significant 
increase in the transitional Reynolds number is obtained by appropriately 
prescribing the wave length and phase velocity of the wall motion. Previously 
developed asymptotic methods are used in the analysis. 



1, Introduction 


Nearly one -half of the total drag of high-subsonic speed, long range 
aircraft, is a result of viscous drag caused by turbulent skin friction* 

Thus, any method that substantially reduces this drag will dramatically 
effect future aircraft design, and the general question of energy conservation. 

In a series of pioneering experiments Kramer, see e.g* [1] , demonstrated 
that significant drag reduction can be achieved in the motion of bodies 
under water by coating their surfaces with compliant materials. Kramer's 
most successful experiments occurred for compliant materials that attempted 
to simulate the skins of dolphins. The precise mechanism for the reduction 
of drag by compliant surfaces is not known. Kramer, and other investigators 
conjectured that the drag reduction resulted from a delay in transition from 
laminar to turbulent flow. However, it is presently believed that this 
effect is unimportant at the large Reynolds numbers in Kramer's experiments. 
Subsequent experiments on drag reduction by Russian scientists, see e.g. [2],^ 
substantiated and developed Kramer's ideas. However, experimental results 
on drag reduction by compliant walls for air flows are inconclusive. A 
critical and detailed survey on compliant wall research is presented in [3] . 

In the past decade, fundamental experimental advances have been made 
in understanding how turbulence is produced in a turbulent boundary layer. 

These experiments show that there is an organized sequence of flow events 
that give rise to boundary layer turbulence. The basic feature of turbulence 
producing flows is the formation and subsequent motion of fluid ’’bursts" 
near the wall. The experiments indicate that most of the turbulence 
and hence most of the drag, is produced during the formation of the bursts 
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and their eventual motion away from the wall. A discussion and summary 
of the burst formation process is given in [ 3 ] • 

Bushnell [3,4] has conjectured that compliant walls produce drag 
reduction by interfering and altering the bursting process. This mechanism 
requires the walls to respond with sufficiently large amplitude, and 
sufficiently high frequency. A variety of compliant walls have been 
employed using different materials and different structural configurations. 

They are described in [3,4,5] . However, none of these walls have provided 
the required motions, or resulted in significant drag reductions for air flows. 

In this report, we present a new mechanism to achieve the desired 
wall motion. It employs slightly curved walls subjected to external pressures, 
and the ensuing nonlinear, large amplitude motion of the wall. That is, it 
relies on the dynamics of snap-buckling of curved elastic plates. To 
demonstrate this mechanism, we consider a shallow elastic arch, or equivalently 
a wide curved elastic plate, whose axis is aligned with the flow direction. 

To simplify the problem, we neglect the interaction of the arch with the 
air flow, and assume that a prescribed static pressure is applied to the 
surface of the arch, A further simplification is made by assuming that the 
undeformed shape of the arch, the pressure, the disturbances, and the response, 
vary sinusoidally in space. Then it is possible to solve the resulting dif- 
ferential equations exactly. 

We show in Section 3, that as the pressure increases from zero, the 
arch assumes stable, static equilibrium states with relatively small ampli- 
tudes. They are called unbuckled states. At a critical value of the 
pressure, the unbuckled state becomes unstable because there are no nearby 
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equilibrium states for pressures slightly greater than the critical value. 
Then due to any small disturbance, the shell jumps, or snap buckles, 
resulting in a large anq}litude periodic motion. This motion is obtained 
in Section 5, in terms of the geometrical and material parameters of the 
arch, and the amplitudes of the fluid-imposed disturbances. 

In Section 6, we apply these results and show that for experimentally 
feasible arch geometries, and for certain materials like mylar, the ampli- 
tudes and frequencies of the motion are within the limits required by 
Bushnell*s hypothesis for drag reducing compliant walls. However, the 
theory we have presented is not yet complete, because it was asstuned that 
the fluid and arch motions are decoupled. Hence this theory does not account 
for the coupling or give a specific mechanism for transferring energy from 
the fluid flow to the arch to maintain the desired periodic motion, A more 
refined theory is required, which accounts for the fluid-solid interaction. 
The asymptotic expansion, in terms of rational functions, of the exact 
solution that is obtained in Section 5, suggests a new asymptotic method 
for analyzing jump phenomena like, snap buckling, combustion, and population 
explosions. Standard as 3 miptotic and perturbation methods are not applicable, 
because in jump phenomena the order of magnitude of the solution changes 
during the jump. The new method, which we call the method of rational 
functions has been presented in [6,7] • 

Experiments with boundary layer flows over flat plates suggest that 
the delay of transition and drag reduction can occur by suitably driving 
the flat plate with a traveling wave, see e,g, [2,8,9 1. We have considered 
the mathematically sinq)ler problem of Poiseuille channel flow with driven 
walls [10 ], The results are summarized in Section 7, where we show that 
transition from the laminar Poiseuille flow can be delayed by suitably 
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varying the phase velocity of the wall. Previously developed [ll,12l 
asymptotic methods were employed in the analysis. 
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2, Formulation of the Dynamic Snapping of an Elastic Arch 


A shallow elastic arch of length L, cross-sectional area A, and cross- 
sectional moment of inertia I, is subjected to a lateral force per unit 
length P(X,T). Here X is the axial coordinate of the arch, 0 ^ X L, 
and T is the time. The vertical position of the undeformed arch is 
Y = (X) , where is a specified function. If the deformed position 

of the arch is denoted by Y = y(X,T), then the "excess" vertical dis- 
placement W(X,T) is defined by 

W(X,T) H w°(X) - y(X,T) . (2.1) 


We now define the following dimensionless variables and parameters. 


X = X/L , t E (E/p)^^^T , w(x,t) = W(X,T)/L 

2 

p(x,t) = ^ P(X,T) , w°(x) = W°(X)/L , 

TT ^ 


( 2 . 2 ) 


where P and E are the mass density per unit length and the Young's modulus 
of the arch, respectively. In terms of these variables, the initial boundary 
value problem for the arch is [13,14], 


w - sw 

XXXX XX 

+ [7t^/(4£)]w^^ = - sw^ + [tt^/(4£)]p , 

(2.3a) 


,1 


S(t) H 

1 (w - 2w^ w )dx , 

X X x'^ ’ 

(2.3b) 

V - w - 0 

XX 

, for X = 0,1 , 

(2.3c) 


w(x,0) = f(x) , w^(x,0) = g(x) , 

where the dimensionless parameter ^ is defined by 


(2.3d) 
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i = lVi 


(2.4) 


The quantity s(t) that is defined in (2«3b) is the dimensionless axial 
stress in the arch. It is independent of x. The boundary conditions (2,3c) 
imply that the ends of the arch, x = 0,1 are simply supported. The prescribed 
functions f(x) and g(x) are the dimensionless initial vertical displacement 
and initial vertical velocity of the arch, respectively. 

We analyze (2,3) for sinusoidal arches, with sinusoidal initial data, 
and static but spatially sinusoidal pressures. Thus we have 


Wq “ h sin TTx , p - q sin ttx , 
f = f^ sin irx , g = Sq sin ttx , 


(2.5a) 

(2.5b) 


where h is the dimensionless maximum arch rise and q is the dimensionless 
maximum pressure. Then (2,3) and (2,5) has solutions 


w(x,t) = A(t)sin irx , 


( 2 . 6 ) 


where A(t) is the excess vertical displacement at the center of the arch, 
and it satisfies the initial value problem, 


A" + aA - 3hA^ + A^ - q = 0 
A(0) - fjj , A' (0) = gp . 


(2.7) 


A prime denotes differentiation with respect to t, and the parameter a is 
defined by 


a = 2(h^ + 21 ) 


( 2 . 8 ) 



3, *rhe Static States 


It follows from (2.7) that the amplitudes of the single mode equilibrium 
positions of the arch satisfy the cubic equation, 

q = aA - 3hA^ + . (3.1) 

We refer to, (3.1) as the static response equation, and to its graph as 
the static response curve. 

The number of real roots of (3.1) depends on the sign of the discriminant, 

D = h^ - a/3 = . (3.2) 

We consider only arches with D > 0. Then there is a unique static state 
for q in the intervals q < q^, and q > q^, see Fig. 1 . However, for q in 
the interval q^ < q < q^» there are three static states. The lower branch 
of this curve, with amplitudes in the range, 0 A f.A^, is called the 
unbuckled branch. The upper branch (A L called the buckled branch. 

The remaining states with amplitudes in the interval, A ^A„, lie on 

the intermediate branch. 

As q increases from zero, the equilibrium states of the arch lie on the 
unbuckled branch. If q exceeds q^, then the arch must jump dynamically 
towards the buckled state. Similarly if the arch is in an equilibrium 
state corresponding to a point on the buckled state, and q is slowly 
decreased below q^^, the arch will suddenly jump towards the unbuckled 
state. These d 3 mamic transitions are called snap buckling. They are 
described by the initial value problem (2,7) for disturbances given by (2.5b), 

Consequently the dimensionless pressure q^ and q^ are called the upper 
and lower buckling loads, respectively. Their values are easily obtained 
from an elementary analysis of (3,1) as. 
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q - 4h2, - 2D‘ 


(3.3) 


= 4hil + , 


3/2 


The corresponding values of the amplitudes are, 

= h - = h + . (3.4) 

The amplitude of the buckled state, -4jj(q^)» corresponding to q = q^ is 
obtained from (3.1) as 


Aj^(q^) = 3h - 2Aj = h + . (3.5) 

The stability of the equilibrium states to single mode disturbances 
of the form (2,5b) is determined by linearizing (2,7) about any equilibrium 
state A = A^. From an analysis of this linearized problem, we conclude 
that the unbuckled and buckled branches are stable, and the intemediate 
branch is unstable. In addition, the limit points, q = q^, A = A^, and 
q = q^ , A = A^ , are unstable equilibrium states. We wish to study the 
motion that results (snapping) when the arch is given a small disturbance 
from the equilibrium state corresponding to a limit point. 
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4, The Phase Plane Analysis 


The qualitative features of the dynamic motion of the arch will now 
be obtained from a phase analysis of (2«7). The system (2,7) describes 
the oscillations of a conservative mechanical system. The phase plane 
diagram is sketched in Figure 2, for D > 0, and q in the interval q^< q < q^» 
where Aj^(q), A 2 (q), and A^(q) are the amplitudes of the equilibrium states 
corresponding to the unbuckled, intermediate, and buckled branches, 
respectively. Each trajectory is characterized by a value of the initial 
energy of the system. 

The equilibrium points A^ and A^ are centers in the phase plane, and 
the unstable equilibrium point A^ is a saddle. The closed curves surrounding 
and lying inside the separatrix in Figure 2 represent periodic 
oscillations of the arch about the unbuckled (buckled) equilibrium state. 

The closed curves lying outside of the separatrix and enclosing all three 
equilibrium points in their interior, represent large amplitude periodic 
motions of the arch swaying back and forth between the unbuckled and buckled 
states. It is the latter, large amplitude motion that is of particular 
interest for drag reduction structures. 

As q q^, the equilibrium states A^ and A^ coalesce to form the limit 
point. In the phase plane, the center and the saddle coalesce to form a 
new singular point, which is a cusp. The resulting phase plane diagram 
is sketched in Figure 3, The closed curves surrounding both equilibrium 
points, and lying outside the separatrix, represent large amplitude, 
periodic motions of particular interest because they correspond to the snap 
buckling of the arch. 
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It Is possible to calculate the local periodic motions about and A^ 
in Figure 2, and A^ in Figure 3 by the Poincare-Linstedt [15 ] perturbation 
method. Since our primary interest is in the large amplitude motion of the 
arch, this method is not applicable. In the next section, we determine 
the snap buckling motion by solving (2,7) exactly in terms of elliptic 
functions. Then an expansion of the resulting solution is obtained for 
small initial disturbances from the limit point. This suggests a new 
asymptotic method. for analyzing jump phenomena. This new method, which 
we call the Method of Rational Functions, is described in [6,7], 
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5* The Dynamics of Snap Buckling 


To analyze the snap buckling of the arch, we assume that q = and 
the arch is in equilibrium at the limit point, A = A small initial 

displacement disturbance e is given to the arch, so that, 

A(0) = A^ + £ , A*(0) = 0 . (5,1) 

The ensuing large amplitude periodic motion describes the snap buckling of 
the archo In our analysis, we assume that e < 0, A similar analysis applies 
for £ > 0* Then the motion is periodic, . but it lies on a closed curve 
inside the separatrix in Figure 3o 

It is convenient to transform the dependent variable in (2,7), to 
analyze the motion relative to the limit point. Thus, we define v(t) by 

A(t) = A^ + v(t) . (5.2) 

By inserting (5,1), (5,2), and q = q^ into (2.7), the initial value problem 
is transformed to, 

v" - + v^ = 0 , v(0) = e , v*(0) = 0 , (5.3) 

where the parameter D is defined in (3,2), The phase plane diagram for 
(5,3) is the same as Figure 3 with the origin shifted to the limit point. 

By integrating the first integral of (5.3), we obtain a representation 
of the solution of (5,3) in terms of the elliptic function 

= COS T (5,4) 

and the roots of a quartic equation, see [ 6l for details. For small z the 
roots are obtained by a perturbation method. This leads to the representation 
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v(t;e) = 


. (5.5) 


s'a - /J) 

1 4,(t) 

L \1 - /3)/ J 

e + O(e^) 

[1 + 4)(t)] + (” - /3) 

/l + /J\ 

1 + f U(t) 

. \\-^J J 

E + O(e^) 


1/2 

where the parameter 6 is defined by 6 = 4 d , and 

1/2 


T = 




(-e) + 


• t 


(5.6) 


The function <|)(t) is periodic of period 4K, where K(k) is the complete 
elliptic integral of the first kind with modulus k. The modulus is given by 


■(^) 


+ 0(e) 


(5.7) 


From tables of the complete elliptic integral of the first kind [16] we get, 

/3 ^ 


K(k) ==» Kq + 0(e) , where Kq = K 


(H 


2.7693 


(5.8) 


It follows from (5.5) and (5.6) that the period of the motion is 

1 /9 1 /9 

tQ = 4/2K/(BE)‘^'^ = y^Yf2 • 


(5.9) 


The mathematical structure of the snap buckling of the arch, and perhaps 

other jump phenomena, where the solution changes from small [0(e)] to 

large [0(1)1, is revealed by the rational functions (5.5). Since 

1 + 4^(1^) - 0 only for t = T^^(e) = 2(1 + 2n)K, n == 0,1,2,..., the leading 

term in the denominator, 1 + 4 >(t), does not vanish for t in the interval, 

T < T , Hence, the denominator changes from 0(1) to 0(e)# and then 

n n+i 

back to 0(1) as t varies through t^. Since the numerator is 0(e), the 
solution jumps from 0(e) to 0(1) and back to 0(e) as t varies through 
The jump occurs smoothly, but rapidly. That is, the solution remains 0(e) 
for most of the period, and jumps to 0(1) and returns to 0(e) over a 
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relatively short interval of time. The period of the motion is long, as 
we deduce from (5.9). 

If the Poincare-Linstedt method is applied to (5.3), the leading term 
is given by (5.5), with only the 0(1) term in the denominator. Thus, the 
approximation given by this method is unbounded at n = 0, !,•••, 

and hence it is invalid. 
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6. Application to Drag Reduction Compliant Walls 


For an arch of rectangular cross-section of breadth b and depth d, 
the parameters £ and D that are defined in (2*4) and (3o2), respectively, 
are given by 

n = (1/12) (d/L)^ , D = J [(H/L)^ - (l/3)(d/L)^l . (6.1) 

Here H = Wp(L/2) is the arch rise, ioe« it is the physical height of the 
center of the undeformed arch above the horizontal* 

In physical variables, the pressure and the static displacement 
at the center of the arch corresponding to the limit point q^, are 


== ^ r/d\ 1 /Hy d\^ 

b 4 bL ^ 4 [VLy3^j\Ly 


+ 


9 


( 6 * 2 ) 


W = Lw (1/2) = LA = L r I- - ] * 

u u u L L J 

We have used (2,2), (2*6), (3*3) and (3*4) to obtain (6.2)* Similarly, 
the physical displacement W^ at the center of the static buckled arch for 
q = is obtained from (3.5) and (2.2) as. 




(6*3) 


The maximum displacement W of the arch during the motion is evaluated from 

m 

(5*5) and (3.4) as 

^ L [S + + O(e^) j * (6.4) 

Finally, the period Tq in seconds of the large amplitude motion, with 
q = q^, is obtained by inserting the time scaling in (2*2) into (5*9)* 

This gives. 
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Tq 2 ^0 2.1/4 1/4 Ve; 

TT 7T J D 


4/2K^L . -1/2 . .. 

° (&) + 0 ( 1 ) 




(6.5) 


Typical dimensions (in inches) of panels used in drag reduction experi- 
ments [3l-[5l are, 

.02 1 L ^ .035 , o002 1 H £ .003 


,0001 < d < .01 


(6o6a) 


Then the parameters (H/L), (d/L) , and D lie in the intervals, 


.07 < H/L < .15 


,005 ^ d/L 1 .286 , 


0 < D < .0075 


(6.6b) 


For (H/L) large compared to (d/L), then D, and the quantities in 
(6.2) -(6,5) are given approximately by 

P _4 , .,-3 

6 >/ 3 ’ 




(6.7) 


^ - 1 + 2 r- 

W ^ H , W, «=• H , W (1 + /3)H 

u /3 » b ’ m ^ 

Thus and are approximately 5 and 6,5 times larger, than W^, respectively 
if (H/L) » (d/L), 


If we denote the height of the arch at maximum displacement by H , 


then we obtain from (6.4) 


. 1/2 


H = H +W = 2H + 3LD + 0(e) 
m m 


( 6 . 8 ) 


By substituting the maximum values of H, L and D given in (6.6) into (6.8), 

we obtain for the maximum value of H 

m 
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(6.9) 


H .015 inches 
m 

Significantly smaller values of are obtained for other values of H, L 
and D in the ranges given in (6.6). For example, if H/L is large compared 
to d/L then using (6.7) in (6.8) we obtain 

(2 + /3)H . (6.10) 

Then by using the minimum value of H from (6.6a), we get H. ^ .0075 inches 

m 

from (6.9). It is conjectured [3] that wall motions with relatively large 

amplitudes are necessary for drag reduction. In experiments conducted at 

Langley with other compliant wall structures and mechanisms, it was difficult 

to excite motions with amplitudes as small as .0008. The present analysis 

suggests that amplitudes ten to twenty times larger can be obtained using 

the snap-buckling of curved structures. 

If the maximum amplitudes of compliant walls are too large, they will 

induce drag because of a "roughness” effect. For turbulent boundary layer 

flows, an empirical rule is that the maximum amplitudes should not exceed 

100 v/U , where v is the kinematic viscosity of the fluid and U is the 

free stream velocity of the flow [21] . For the typical experimental values 
“A 2 

of V = 1.5 X 10 ft /sec. and in the range, 3 ft /sec ^ 25 ft /sec, 

the empirical rule requires that the maximum amplitudes lie between .0072 
inches and .06 inches for a hydraulically smooth surface. The geometrical 
parameters H, L and D can be selected so that the maximum amplitudes lie 
within these values. 
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It has been conjectured from experiments and a study of data, see e*g, 

-I 

[3lj that wall frequencies 2000 Hertz are necessary for drag reduction* 

For this frequency, we conclude from (6*5) that the corresponding amplitude 
c of the disturbance is given approximately by, 

-e 582 ^ 10^(p/E)L^/D^^^ * (6.11) 

In Table 1, we give the values of -e obtained from (6.10) for typical 
materials used in drag reduction studies® The values of p/E were obtained 
from [ 5l e 

To interpret the results in Table 1, we recall that c represents a small 

displacement disturbance applied to the arch. We observe from the table 

that at a frequency of 2000 Hertz, mylar would be a suitable material for 

2 1/2 

a wide range of the geometric parameter L /D « That is, small disturbances 

can produce a periodic motion at the required frequency. The materials 100 PPI 

2 1/2 

foam and RTV are suitable provided that L D is sufficiently small. Other 
materials such as PVC that have been used in drag reduction experiments are 
unsuitable because they require a relatively large value of -e to excite 
the periodic motion with the required frequency. 

Criteria for selecting the magnitude of e are not given by the present 
theory. The parameter z crudely represents disturbances produced by air 
flow over the arch, such as turbulent boundary layer fluctuations. If the 
magnitude of these fluctuations is known, or can be measured in an experi- 
ment, then, by using the analysis and results of the present investigation', 
it is possible to crudely estimate the geometric and material properties 
of the arch that are necessary to produce the desired frequency and 

amnlitude of resnonse. These results apply only to Initial disturbances of the 
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form (2.5b), If more complicated disturbances occur then the response will 
contain more than one mode. Additional modes introduce more equilibrium 
states of the arch than are illustrated in Figure 1, see e.g* [14], 

Aerodynamic and structural damping have been neglected in the formu- 
lation (2.3). The presence of a small amount of structural damping will 

change the center at A« in Figure 3 to a stable spiral point. Then instead 

3 

of snap buckling occurring as a periodic oscillation, there will be a 
damped oscillation with the solution approaching the buckled equilibrium 
state at A^ as t , That is, energy is lost by the system due to structural 
damping, and the arch eventually comes to rest at the buckled equilibrium 
point. In order to maintain a periodic oscillation, energy must be supplied 
to the arch from its environment, to replace the energy lost by dissipation. 
This energy may be obtained from the pressure induced by aerodynamic forces, 
and from turbulent boundary layer fluctuations. For example, the turbulent 
boundary layer fluctuations may provide sufficient disturbances to the arch 
at regular intervals to continually excite a periodic, or nearly periodic, 
large amplitude oscillation. However, the aerodynamic forces may also damp 
the motion, A more sophisticated study is required to determine whether the 
net energies supplied to the arch by the aerodynamic pressures and turbulent 
boundary fluctuations are sufficient. Then it is necessary to consider the 
coupled aerodynamic-structural interaction problem. 
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7, Poiseuille Channel Flow with Driven Walls 

The flow of a viscous, incompressible fluid between two horizontal, 
rigid and infinite walls that is driven by a horizontal pressure gradient 
is called Poiseuille channel flowo It is described mathematically by an 
initial -boundary value problem for the two-dimensional Navier-Stokes 
equations. Plane Poiseuille flow is a laminar flow solution of this ppoblem 
for all values of R. The Reynolds number that characterizes the flow is 
R = ’ here, 2H is the channel gap , ^ is the kinematic viscosity and 

Uq is the centerline velocity of this flow. 

It follows from the linearized theory of hydrodynamic stability [17] 
that plane Poiseuille flow is stable (unstable) for R < R^ (> • The 

critical Reynolds number, R^, is determined numerically as, 5772, 

see e.go tisl. In addition, a branch of unstable spatially and temporally 
periodic solutions of the Navier-Stokes problem bifurcates subcritically 
from the laminar flow solution at R = R [19]. Thus, R = R is a bifur- 
cation point of plane Poiseuille flow. 

In [10] we have modified this classical channel flow problem by driving 
the upper wall with a specified vertical motion to determine if transition 
can be delayed to a Reynolds number ^ R^o In the dimensionless variables 

X H x/H , y = Y/H , t = (Uq/H)T , (7.1) 

the wall motion is given by 

y = 1 + 6G(x,t) , G(x,t) = cos(ax-ct) . (7.2) 

Thus, 5G(x,t) is a traveling wave of amplitude 6 and phase velocity c/a, 
where a and c are the dimensionless wave niimber and circular frequency, 
respectively. For simplicity of presentation, we assume that 5^0. In 
* 

This H should not be confused with the maximum arch rise of Sections 2-6. 
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terms of a dimensionless stream function ’l'(x,y,t), the two dimensional 
Navier -Stokes equations are reduced to 

AT _ i + J(T,AT) - 0 . (7.3a) 

u R. 

2 

Here ^ and A are the two dimensional harmonic, and biharmonic operators, 
respectivelyo The nonlinear operator J is defined for any two functions f, 
and g by 


J(f,g) = f g - f g • (7.3b) 

' y X X y 

The boundary conditions on the walls are, 

H'^(x,-l,t) = T^(x,-l,t) = 0 (7.3c) 

'l^^Cx,! + 6G(x,t),t) = -6G^(x,t) , ¥y(x,l + 6G(x,t),t) = 0. (7,3d) 

The boundary conditions as |xj ->■ ■» are, 

’F(x,y,t), and its derivatives, are bounded as |x[ ” • (7.3e) 


For the dependence of the wall motion on the Reynolds number, we assume 


that , 


c 




(7.3f) 


where c^ and w are parameters to be specified. For simplicity of presentation, 
the wave number a is assumed to be independent of R. The perturbed problem 
that we analyze is to determine solutions of (7.3) that are periodic in x 
with wave length and periodic in t with period 2tt/c, The perturbed 

problem is a forced oscillation problem, where the wall motion is the 
"forcing function". Thus, we are seeking solutions with the same period 
as the forcing function. 
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The classical hydrodynamic stability problem of Poiseuille channel 
flow is obtained by setting 5 = 0 in (7*3) • We refer to the problem of 
determining periodic solutions of periods 2 ' u/cl in x, and 27r/(i0 in t, of 
(7.3) with <S = 0 as the bifurcation problem* Here, the wave number a, and 
the circular frequency to are to be determined* The bifurcation problem 
possesses the laminar flow solution, 

¥ = E y . y3/3 , (7.4) 


for all values of R, and a = (o = 0* It is the classical plane Poiseuille 
flow* In the terminology of bifurcation theory, it is the basic solution of 


the bifurcation problem* 

The stability of the basic solution (7*4) is determined by linearizing 
(7.3) with 5=0, about the basic solution. Solutions of the resulting 
problem are obtained in the form. 


-Ct ikx 

'{' = e ® 4>(y) 


(7.5) 


where the wave number k is real* Then C, and cf) are determined as the eigen- 
values ?j(k,R), and the eigenfunctions, <j>j(y;k,R) of the Orr-Sommerfeld 
problem [17]. 

We have employed a modified version of the Chebyschev method [18] to 
solve the Orr-Sommerfeld problem numerically. The resulting critical Reynolds 
number R^, the wave number k = k^, and the frequency co “ criticality 

are estimated as. 


R 5772*22 , k 1*02056 , w e im .264002 * (7.6) 

c ’c cl 

Here i^j^(k,R) and ‘l^ 2 ^(y;k,R) are the eigenvalue and eigenfunction corresponding 

to the lowest mode, and Re dimensional variables, the 

wavelength at criticality is approximately three times the channel depth* 

Thus, the laminar flow is stable for R < R , and it is unstable for R > R . 

c’ c 
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A generalization of the method presented in [11,12] was employed [10] 
to determine asymptotic expansions of the solutions of the perturbed problem 
(7o3) as 5 0, These expansions are obtained as perturbations of the solu- 

tions of the bifurcation problem. Since these perturbations are singular, 
the method of matched asymptotic expansions is used in the analysis. 

We now give a brief summary of the results of this analysis. All details 
of the analysis are omitted; see [10] for details. 

On the results for the resonance case, which is defined by. 


a = k 


, c(R^) = Cq = 


(7.7) 


are discussed in this report. The resulting solutions of the perturbed 
problem in the non-resonance cases are in general, quasi -periodic functions; 
see [10] for details. Asymptotic expansions of the solutions of the perturbed 


problem for R near and in the resonance case are given by. 


4^ = + 2r(R)Re|e 


i [k x-(o) -hw^6^^^)t-0+O( 6) t ] 


.1/3 


. 2 / 3 , 


4>i(y)^6 + 0(6 , 

(7.8) 


where and are given in (7.6), w is defined in (7.3f), 4>^(y) is the 
Or r -Sommer f eld eigenfunction, and ? Is defined by 


R - + 0(6) . (7.9) 

Thus, by varying K and w, the Reynolds number and the forcing frequency c 
are varied respectively near their critical values. The amplitude r, and the 
phase shift 9 in (7.8) are the real roots of the amplitude equations. 


+ p5r = ±(Y^ - , (7.10a) 

0 = 0 + . (7.10b) 
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The quantities p, and q in (7,10) are linear functions of w that are numeri- 
cally evaluated from integrals involving the Or r -Sommer f eld eigenfunction 
and its adjoint eigenfunction. They are given by, 

1.466825 - (,0056057)w , q^ -.582325 + (,001004)w . (7.11) 

The constants Y and a in (7.10), which are evaluated similarly, are given by 

Yc?? .000163 , a ^-.229238 . (7.12) 

The response (7.8) is a small perturbation of the laminar flow 

1/3 

The amplitude of the perturbation is proportional to r5 • Thus, the small 
disturbance, whose amplitude is proportional to 6, produces a much larger 
response. The exponential factor in (7,8) gives the propagation of the 
transverse normal mode down the channel. Thus, the perturbation prop- 
agates with the same wave length in x and period in t as the impressed wall 
motion, but it is phase shifted by 6(R), see (7,3f) and (7.8). We refer 
to the solutions (7.8) as perturbed laminar states. 

An elementary analysis of the amplitude equations (7,10) which we do 
not present, shows that there are four critical values of w. They are 

denoted by w_ , w , w„, and w . The critical values w and w are the 
L p H p 

zeros of p = 0 and q = 0, respectively. The critical values w < w are 

L H 

the two roots of 

p^(w) - 3q^(w) = 0 . (7.13) 

By using the values in (7.11) and (7.12), these critical values are given 
approximately by, 

Wj^<5« 116 , 262 , Wgtrs? 338 , w*«? 585, • (7.14) 

The critical values define five intervals of values of w where the response 
curves have different qualitative features. Typical response curves 
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for each of these intervals are sketched in Fig. 5. The perturbed laminar 

states have a maximum amplitude at R = R^, The value of this maximum 

and its location depend on w, as illustrated in Fig. 4. Asw->w,r ^ „ 

m 

* 

and R^ ->•«•. Thus, when w = w , the perturbed states give a continuous 

transition from the laminar state to the bifurcation state, as shown in 

Fig. 4e, The dotted curve in each of these figures is the solution of the 

bifurcation problem that branches from R = R^, 

The significance of the critical value w is that R > R (c R ) if 
® p m c c'^ 

w < w_ (>w). Atw=w,R and R coincide. For w in the open interval 
p ' p' p’ m c ^ 

(w w ), the response curve gives the amplitude as a single valued function 

3 £1 

of Ro For w outside this interval the response curves are multiple-valued. 

At w = w , and at w = w , the response curve has a vertical slope near 
L £1 

R = R^. Multi-valued response diagrams define two additional critical 
values of R« They are denoted by R^, and R^ as shown in Fig. 4a. They are 
characterized as the points of vertical slope. Their values, which depend 
on w, are obtained directly from (7,10a). 

The interpretation of the results depends on the stability of the 
solutions. For many nonlinear stability problems it can be shown that 
solutions are stable (unstable) if 

>0 (<0) . (7.15) 

That is, they are stable (unstable) if the amplitude of the response increases 
(decreases) with the Reynolds number. We shall discuss the results assuming 
that this statement is true. However, it is important to emphasize that the 
present analysis determines periodic solutions near the laminar state that 
have the same wavelength and frequency as the wall motion. They may be 
other periodic solutions of the perturbed problem, such as subharmonic 
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solutions o In addition, there may be other periodic, and non-periodic 
solutions of the initial value problem for the system (7«3). Since the 
precise long-time behavior of the solutions of this system is unknown, it 
is difficult to assess the significance of these additional solutions. In 
particular, the stability of a class of solutions of (7,3) was investigated 
in [20 ]o There the frequency of the wall motion was fixed at criticality, 
so that in the present notation, only the special case w e 0 was considered. 
The Stuart -Wat son method was employed leading to a first order, ordinary 
differential equation for the slow time variation of the complex valued 
amplitude. This equation possesses time independent solutions, which were 
shown to be stable for R < R- < R , and unstable for R > R- , Thus, if 
the Stuart-Watson method provides a valid approximation to solutions of 
(7,3), then the analysis in [20] suggests the possibility of the bifurcation 
of additional periodic solutions at R = R^, An investigation of the stability 
of the solutions (7,8) is currently in progress. 

We first discuss the results for w < w ; see Figs, 5a, b. Then as R 

P 

increases from zero, the fluid is in a stable, perturbed laminar state. The 
amplitude of the wave increases until R = R^, Any further increase in R 
results in an unstable state. This suggests that as R increases through R^ 
there is dynamic transition to another state of the fluid. This transition 
may be rapid and involve large amplitudes. However, if w < w^^, and R = R^, 
then the solutions may execute a "small" amplitude jump from the upper branch 
in Fig, 5a to the stable intermediate branch with r = r^. Then by further 
increasing R, a second transition will occur at R = R^ to some other state 
of the fluid (possibly turbulent). It is important to boserve that for 

w < w , R , and hence R , exceed R , So that the jumps we have described 

p m u c -1 r- 

occur when R exceeds R^, That is, transition is delayed. The maximum 
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excursion of to the right, and hence the largest delay occurs when 

w 180, as is deduced from our analysis. The corresponding value of 

is denoted by R , and it is given by 
m 

^c*-^ 28.0636^'^^ + 0(6)1 . (7.16) 

Then for the typical experimental value of 6 = .002, we obtain from 

(7.16) that R 8344. This gives a 45% increase in the transitional Reynolds 
m 

number. Thus, wall motions provide a method for nonlinear stabilization. 

This interpretation may be of practical significance for the control of 

turbulence in channel, pipe and other viscous flows. Additional control 

of the transition might be achieved by oscillating the lower wall so that 

it is phase shifted with respect to the upper wall. 

The stability of the fluid in the presence of the wall motion is 

determined by two competing physical effects. The wall motion disturbs the 

laminar flow, which tends to destabilize the fluid. However, the resulting 

motion is out-of-phase with the wall motion. If the phase shift is small 

then the disturbance effect of the wall wins, and there is early transition. 

However, if the phase shift is sufficiently large, then the instability is 

restrained by the wall motion. This is confirmed by the variation with w 

of the phase shift corresponding to R == R , see Fig, 5, The phase shifts 

m 

are larger when R exceeds R , i.e. when w < w , 

^ m c’ p 

A paper describing the results of this investigation has been submitted 
for publication. 
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8. Concluding Remarks 


Two mechanisms for drag reduction _^in the flow over a flat plate art. 
considered. The first mechanism employs compliant walls and Bushnell’s 
conjecture that compliant walls produce drag reduction by interfering 
with and altering the formation and motion of the fluid bursts in the 
turbulent boundary layer adjacent to the wall. Experiments indicate that 
most of the turbulence, and hence most of the drag associated with a 
turbulent boundary layer, is produced during the formation of these bursts 
and their motion away from the wall. The walls must respond with suf- 
ficiently large amplitude and sufficiently large frequencies to produce 
the required interference. The proposed mechanism employs slightly curved 
walls and relies on the large amplitude motions due to snap— buckling of 
curved elastic structures. The mechanism is demonstrated for a simpler 
problem for a shallow, curved, elastic arch. For a suitable pressure dis- 
tribution the dynamic problem of the arch is solved exactly. Employing 
this exact solution, we show that for experimentally feasible arch geometries 
and for wall materials like mylar, the amplitude and frequencies of the 
motion are within the limits required by Bushnell’s hypothesis. Our 
analysis of the snap-buckling mechanism is incomplete since the fluid and 
structural motions are assumed to be decoupled. A more refined theory 
that accounts for the interaction of the fluid with the structure is 
required to describe a mechanism of transferring energy from the fluid flow 
to the structure to maintain the desired wall motion. 

The second mechanism is concerned with delaying the transition from 
laminar to turbulent flow along a wall by driving the wall. We have 
considered the related but simpler problem of Poiseuille channel flow 
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whose walls are driven by a periodic traveling wave. It is shown that 
by appropriately selecting the wave number and phase velocity of this 
traveling wave a significant delay in the transitional Reynolds number 
can be achieved. This delay occurs because of the phase shift of the 
resulting periodic fluid motion with respect to the prescribed wall motion. 
Additional investigations are required to determine if a similar delay 
in transition can be achieved for boundary layer flows over a driven wall. 
Recent experiments [8] suggest that this may be possible. 
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TABLE 



*100 PPI Foam is a polyurethane foam. 
**RTV is a silicone rubber. 


-e X 10^ inches 


= .01 

l2/d1/2 ^ Q5 

= .1 

.98 

4.9 

9.8 

2,088 

10,440 

20,880 

8,350 

41,750 

83,500 
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